Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  962 1776 (N38E31 L= 98)  2726    506        8 ( 841/1651-1046/1872) -23.7 E46-E21 N29-N44
      2  472 1672 (N18E59 L= 70)  3558    525        9 ( 383/1617- 561/1765)  33.7 E66-E50 N13-N25
      3  541 1539 (N11E45 L= 84)  1637    362        7 ( 459/1458- 615/1599)   8.1 E54-E40 N06-N14
      4  246 1507 (N04E75 L= 54)  1222    207        2 ( 220/1474- 271/1535)  45.0 E82-E70 N01-N05
      5  478 1463 (N05E47 L= 83)   298     73        1 ( 456/1442- 505/1480)  45.0 E49-E44 N03-N06
      6 1338 1329 (N20W09 L=139)   617    123        2 (1288/1295-1390/1366) -15.9 W06-W13 N16-N22
      7 1562  539 (S16W49 L=179)   907    189        3 (1475/ 487-1644/ 588)  57.5 W45-W54 S21-S10
      8 1568  440 (S20W56 L=185)   530    105        2 (1535/ 407-1595/ 468) -45.0 W52-W58 S22-S18
      9 1133  301 (S44W37 L=166)  1380    281        3 (1012/ 270-1250/ 343)  22.4 W27-W44 S47-S40
  

    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