Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1614 1445 (N21W49 L=194)   720    186        1 (1582/1368-1645/1528)  68.2 W45-W51 N16-N26
      2 1522 1391 (N20W40 L=185)   629    187        0 (1480/1328-1589/1418)  11.3 W35-W45 N16-N21
      3  995 1162 (N15W00 L=145)   388    105        0 ( 973/1108-1019/1216) -90.0 W00-W01 N11-N18
      4 1472  976 (S04W27 L=172)   502    121        1 (1415/ 958-1529/1001) -29.7 W23-W30 S05-S01
      5 1834  906 (S16W56 L=201)   574     90        2 (1804/ 874-1854/ 941)  18.4 W52-W58 S17-S14
      6 1512  731 (S20W27 L=172)  1332    369        4 (1416/ 559-1572/ 904)  55.6 W19-W32 S29-S09
      7 1822  834 (S20W54 L=200)   663     84        2 (1799/ 808-1844/ 865) -51.3 W52-W57 S22-S17
      8 1133  601 (S21W00 L=145)   585    227        0 (1110/ 498-1147/ 700)  74.1 E03-W02 S28-S14
      9  646  481 (S23E35 L=110)   929    232        0 ( 571/ 434- 711/ 537)  19.7 E42-E28 S25-S19
  

    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