Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1218 1571 (N35W28 L=171)   379     89        1 (1179/1544-1270/1589) -18.4 W25-W31 N33-N36
      2 1016 1532 (N38W12 L=155)   222     64        0 ( 998/1505-1032/1562)  36.9 W10-W14 N36-N39
      3  572 1382 (N38E24 L=119)   867    148        2 ( 520/1323- 631/1430) -41.6 E28-E19 N33-N41
      4  334 1161 (N27E45 L= 98)  1638    709        4 ( 283/ 954- 403/1314) -70.8 E50-E41 N14-N37
      5 1781 1171 (S03W54 L=197)   583    115        2 (1743/1136-1819/1196) -39.8 W51-W57 S05-N00
      6  858 1152 (N18E07 L=135)   504    175        1 ( 777/1117- 953/1193) -35.0 E11-E01 N14-N21
      7 1101  929 (N00W03 L=145)  1017    249        3 (1019/ 815-1156/1029) -90.0 W00-W05 S06-N07
      8 1714  935 (S14W43 L=186)   727    105        2 (1696/ 871-1738/ 988)  63.4 W40-W45 S18-S11
      9 1067  639 (S17E06 L=137)   749    138        0 ( 996/ 583-1107/ 675) -54.5 E09-E03 S20-S13
     10  499  349 (S26E56 L= 87)   452    130        2 ( 483/ 312- 516/ 388)  54.5 E58-E53 S29-S22
     11  853  314 (S35E30 L=113)  1091    152        2 ( 810/ 260- 890/ 369)  32.5 E36-E25 S38-S30
  

    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