Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1578 1556 (N25W54 L= 78)   831    121        1 (1519/1530-1624/1595) -45.0 W51-W59 N22-N28
      2  515 1388 (N38E30 L=354)  1281    281        2 ( 426/1318- 660/1462) -31.8 E39-E18 N31-N44
      3 1131 1284 (N20W12 L= 36)   266    131        0 (1106/1235-1180/1320)  31.0 W09-W15 N17-N21
      4  194 1115 (N25E63 L=321)   808    106        0 ( 163/1073- 223/1162) -29.1 E68-E58 N22-N27
      5 1647 1124 (S01W41 L= 65)   281    100        0 (1598/1103-1685/1162) -33.7 W38-W44 S02-N02
      6 1140  819 (S07W03 L= 27)   291    140        0 (1098/ 779-1201/ 838)   0.0 W00-W06 S08-S05
      7 1827  777 (S27W55 L= 79)   429     62        1 (1803/ 748-1850/ 803)  33.7 W52-W58 S28-S25
      8  443  613 (S08E47 L=337)  1440    363        5 ( 371/ 494- 500/ 778)  48.8 E57-E38 S15-N01
      9  281  579 (S09E66 L=318)  2230    240        4 ( 222/ 492- 334/ 662) -80.5 E71-E60 S15-S03
     10  311  621 (S06E59 L=325)   920    125        4 ( 273/ 580- 355/ 657)   5.7 E63-E53 S09-S04
     11  779  507 (S20E26 L=357)   387    159        1 ( 745/ 463- 830/ 565) -76.0 E27-E24 S23-S15
     12 1455  429 (S42W19 L= 43)  1068    107        2 (1422/ 370-1476/ 495)  78.7 W15-W20 S46-S36
     13  474  442 (S20E53 L=331)   799    171        2 ( 425/ 414- 553/ 475) -15.3 E57-E46 S22-S17
     14  661  352 (S29E42 L=342)   741    146        3 ( 635/ 281- 675/ 428)  45.0 E48-E37 S34-S23
     15  961  346 (S35E21 L=  3)   446     88        2 ( 919/ 320- 993/ 384)  23.2 E24-E17 S37-S32
     16 1108  339 (S39E11 L= 13)   410     68        0 (1073/ 308-1133/ 370)  26.6 E13-E07 S41-S36
  

    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