Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1458 1499 (N28W40 L= 29)   690    189        0 (1426/1432-1510/1581)  84.3 W38-W42 N23-N33
      2 1624 1539 (N27W57 L= 47)   509     82        0 (1584/1513-1659/1559) -33.7 W53-W60 N24-N28
      3 1780 1477 (N19W72 L= 61)  1274    125        0 (1744/1422-1824/1533) -76.0 W69-W75 N15-N23
      4  700 1388 (N33E18 L=332)  1032    212        1 ( 656/1286- 749/1473) -66.8 E21-E14 N25-N39
      5 1065 1382 (N28W08 L=358)   247    107        0 (1030/1347-1106/1418) -56.3 W06-W10 N24-N30
      6 1682 1148 (N03W46 L= 35)   807    106        3 (1632/1110-1725/1196) -59.0 W41-W48 N00-N05
      7 1813 1070 (S05W57 L= 47)   502     69        0 (1795/1027-1826/1106) -90.0 W55-W58 S07-S02
      8  362  689 (S07E51 L=298)  1098    196        1 ( 316/ 603- 397/ 769)  45.0 E58-E46 S12-S02
      9 1105  653 (S18E00 L=349)   278    136        0 (1040/ 630-1156/ 685)   8.1 E04-W03 S18-S16
     10  890  428 (S30E19 L=330)   272     93        0 ( 879/ 380- 905/ 476)  76.0 E20-E17 S34-S26
  

    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