Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1080 1615 (N43W12 L=131)   259     40        0 (1053/1598-1107/1629)   0.0 W10-W15 N42-N44
      2 1512 1592 (N35W48 L=166)   763    133        2 (1486/1540-1564/1628)  13.0 W42-W55 N32-N37
      3 1652 1312 (N15W48 L=166)  1245    264        6 (1593/1219-1705/1416) -90.0 W43-W52 N09-N22
      4  227 1278 (N26E64 L= 54)  1203    214        3 ( 169/1251- 317/1320)  -7.1 E75-E51 N24-N28
      5 1370 1241 (N15W25 L=143)   236     97        0 (1345/1203-1398/1280) -78.7 W23-W26 N12-N17
      6  159 1205 (N21E73 L= 45)  1307    102        4 ( 136/1151- 186/1251)  90.0 E78-E68 N18-N24
      7  847 1110 (N12E10 L=108)   227     45        2 ( 820/1088- 877/1131)   0.0 E11-E08 N11-N13
      8 1753  863 (S13W50 L=168)   430     77        0 (1730/ 838-1778/ 885) -53.1 W47-W52 S15-S11
      9  784  378 (S35E26 L= 92)   862    203        1 ( 734/ 315- 839/ 480) -90.0 E28-E23 S41-S27
     10 1205  455 (S34W07 L=125)   350    103        1 (1150/ 442-1255/ 470) -14.0 W02-W10 S35-S32
     11 1136  206 (S57E02 L=116)   453     72        0 (1106/ 193-1175/ 219)   8.1 E05-W02 S58-S54
  

    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