Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1212 1849 (N48W52 L=136)   569     51        1 (1182/1821-1243/1870)   5.2 W47-W58 N46-N50
      2 1100 1820 (N51W40 L=123)   496     55        1 (1079/1789-1123/1847)  12.5 W36-W45 N48-N52
      3  427 1526 (N51E35 L= 48)   367     68        2 ( 408/1495- 450/1554) -51.3 E37-E32 N48-N53
      4  715 1444 (N37E09 L= 75)   349     98        2 ( 671/1417- 761/1467) -39.8 E11-E05 N33-N38
      5  367 1426 (N45E41 L= 42)  1020    182        3 ( 286/1387- 435/1446) -14.0 E52-E32 N40-N47
      6  324 1143 (N27E44 L= 39)  1190    264        2 ( 281/1073- 411/1239) -56.3 E48-E37 N21-N33
      7  355 1137 (N26E42 L= 41)   274     49        2 ( 340/1107- 372/1170)  63.4 E43-E39 N24-N28
      8 1349  951 (S08W16 L=100)   447    115        1 (1329/ 884-1379/1003)  50.2 W13-W19 S11-S05
      9  421  962 (N15E39 L= 44)   471     65        1 ( 402/ 920- 438/1001)  78.7 E40-E36 N12-N17
     10 1321  679 (S23W09 L= 92)   552    157        0 (1305/ 612-1354/ 755)  77.5 W06-W10 S27-S18
     11 1866  707 (S39W61 L=144)   689     68        1 (1833/ 692-1898/ 736)  -7.6 W54-W69 S41-S37
     12 1561  620 (S34W25 L=108)   317     50        0 (1543/ 590-1579/ 648) -78.7 W23-W26 S36-S31
     13 1113  554 (S25E07 L= 76)   207     63        0 (1072/ 542-1145/ 564) -14.0 E09-E05 S26-S23
     14  526  306 (S27E61 L= 22)  2317    246        3 ( 472/ 250- 600/ 367)  10.6 E80-E48 S30-S23
     15 1327  289 (S52E02 L= 81)   692     67        1 (1281/ 271-1366/ 311) -39.8 E05-W01 S55-S49
     16 1222  196 (S58E17 L= 66)   673     35        2 (1194/ 173-1246/ 218) -80.5 E20-E13 S60-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