Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1237 1622 (N36E05 L= 82)  1795    482        3 (1164/1414-1339/1818) -46.1 E18-W07 N25-N51
      2 1011 1545 (N24E16 L= 72)  1445    483        2 ( 858/1388-1172/1642) -10.9 E27-E01 N18-N28
      3 1576 1296 (N26W26 L=114)   809    218        1 (1504/1248-1619/1376)  68.2 W20-W29 N22-N32
      4 1849 1311 (N38W55 L=143)  1002    118        1 (1817/1260-1894/1341)  29.4 W49-W65 N32-N41
      5  698 1226 (S02E24 L= 64)   270     81        0 ( 665/1201- 729/1249)  45.0 E25-E22 S04-N00
      6  165 1204 (S15E68 L= 20)   705     74        2 ( 142/1171- 194/1236) -23.2 E72-E62 S16-S13
      7 1501 1222 (N19W22 L=110)   288     69        0 (1459/1201-1530/1235)  11.3 W19-W24 N17-N20
      8  777 1178 (S03E18 L= 70)   247     27        1 ( 756/1153- 793/1206) -56.3 E18-E16 S04-S01
      9  755 1035 (S12E15 L= 72)   419    147        2 ( 681/ 993- 796/1082)  60.3 E18-E13 S15-S08
     10  393 1034 (S20E40 L= 48)   236     76        0 ( 365/1016- 434/1052)  -9.5 E42-E36 S21-S19
     11 1415  896 (S02W25 L=113)   913    230        4 (1381/ 798-1452/1010) -45.0 W21-W30 S05-N04
     12 1722  777 (N02W51 L=138)   970    312        2 (1671/ 645-1748/ 887) -40.2 W44-W58 S04-N07
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y