Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  964 1547 (N30E04 L=  3)   766    196        0 ( 881/1491-1048/1606) -36.0 E10-W01 N25-N33
      2  246 1332 (N17E56 L=311)   702    100        0 ( 231/1290- 265/1383) -80.5 E58-E54 N14-N20
      3 1664 1341 (N17W44 L= 51)   400     86        1 (1626/1322-1699/1358)  -9.5 W40-W46 N15-N17
      4  147 1205 (N10E67 L=300)   713    121        0 ( 112/1186- 173/1221) -12.5 E72-E62 N08-N10
      5 1766 1136 (N04W51 L= 58)   350     56        1 (1738/1118-1792/1158)  11.3 W48-W53 N03-N05
      6 1854  954 (S06W60 L= 67)   545    104        0 (1837/ 905-1882/1001)  45.0 W58-W63 S08-S03
      7  619  887 (S11E25 L=342)   732    198        0 ( 588/ 780- 650/ 983)  67.4 E28-E22 S17-S05
      8 1045  851 (S14W01 L=  8)   236     59        0 (1006/ 833-1077/ 869) -14.0 E01-W03 S14-S12
      9  987  667 (S25E03 L=  5)   239     40        0 ( 962/ 650-1014/ 689) -14.0 E04-W00 S26-S23
     10  390  616 (S27E46 L=321)   573     83        1 ( 370/ 582- 423/ 646)  31.0 E48-E43 S29-S25
  

    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