Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  900 1461 (N15E18 L=238)  1223    304        1 ( 851/1335- 943/1601) -51.3 E25-E12 N07-N22
      2  338 1442 (N03E57 L=199)  1364    333        1 ( 235/1403- 430/1480)   9.5 E66-E47 N00-N05
      3 1506 1321 (N22W21 L=277)   344     80        0 (1488/1275-1518/1363) -71.6 W19-W22 N19-N25
      4 1041 1255 (N06E05 L=251)   601    149        1 ( 965/1205-1102/1311) -16.7 E10-W00 N04-N07
      5  216 1107 (S19E57 L=199)  2249    466       11 ( 160/1028- 266/1220) -51.3 E64-E50 S23-S12
      6 1508 1041 (N06W27 L=283)  1292    261        0 (1467/ 915-1549/1174) -47.5 W22-W33 N00-N12
      7  661  860 (S25E18 L=238)   821    281        1 ( 532/ 826- 784/ 911)  33.7 E26-E11 S29-S19
      8  265  780 (S38E49 L=207)  1636    491        2 ( 169/ 695- 393/ 817)  25.6 E60-E37 S44-S33
      9 1648  408 (S22W64 L=320)  1808    203        2 (1580/ 365-1695/ 473)  60.9 W59-W68 S26-S17
  

    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