Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  534 1713 (N45E51 L=192)   506     99        1 ( 506/1689- 565/1737)  -7.1 E55-E47 N42-N46
      2 1592 1585 (N42W45 L=288)   589    157        1 (1545/1561-1639/1605)  11.3 W39-W50 N40-N44
      3 1603 1174 (N14W38 L=281)   309     72        1 (1575/1153-1633/1212) -45.0 W36-W40 N13-N16
      4 1785 1091 (N10W54 L=297)   864    161        1 (1757/1014-1808/1162) -71.6 W51-W57 N05-N14
      5 1881 1028 (N06W67 L=310)   993    234        2 (1847/ 974-1931/1062)  20.0 W61-W75 N02-N08
      6  983  887 (S06E02 L=241)   837    202        1 ( 937/ 779-1007/ 985) -71.6 E05-W00 S12-N00
      7 1584  655 (S18W41 L=284)   377     74        0 (1549/ 631-1608/ 683)  45.0 W38-W43 S19-S16
      8 1401  598 (S23W28 L=271)   246     71        0 (1369/ 580-1436/ 622)  -9.5 W25-W31 S23-S21
      9  359  533 (S35E50 L=193)  1824    453        5 ( 254/ 508- 443/ 571)  12.3 E64-E41 S38-S31
  

    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