Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  633 1681 (N27E41 L=195)   800    147        1 ( 603/1619- 664/1762) -37.6 E49-E35 N23-N33
      2  329 1507 (N12E58 L=177)   510     73        1 ( 304/1484- 355/1538) -20.6 E63-E55 N10-N14
      3 1241 1368 (N17W07 L=243)   820    200        3 (1143/1335-1339/1389)   8.7 W00-W13 N14-N19
      4  296 1120 (S11E50 L=186)  1697    458        3 ( 165/ 952- 422/1253)  42.3 E63-E41 S22-S02
      5  184 1120 (S12E62 L=174)   909    154        2 ( 134/1051- 210/1172) -32.5 E69-E58 S16-S09
      6  323 1125 (S11E47 L=188)   402     97        0 ( 288/1077- 352/1160)  51.3 E50-E45 S13-S08
      7  697 1105 (S07E20 L=215)   341     99        1 ( 661/1059- 724/1155) -51.3 E23-E18 S09-S04
      8  142 1084 (S15E67 L=169)  1314    131        5 ( 105/1041- 186/1140) -40.6 E73-E60 S17-S11
      9 1649  906 (S01W41 L=277)   228     32        1 (1622/ 890-1672/ 925)   0.0 W38-W42 S01-N00
     10  221  738 (S35E60 L=176)  1245    103        1 ( 186/ 699- 266/ 790)  40.6 E66-E53 S37-S31
     11 1493  336 (S37W57 L=293)  1445    214        3 (1428/ 303-1580/ 397)  13.2 W47-W66 S39-S33
     12 1244  289 (S47W40 L=276)   541     81        1 (1206/ 268-1291/ 313)  -5.7 W35-W45 S48-S45
  

    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