Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1611 1721 (N28W75 L= 44)  1294     99        2 (1569/1692-1649/1756) -80.5 W70-W81 N25-N31
      2 1536 1672 (N29W60 L= 28)   660     47        1 (1512/1644-1564/1700)  31.0 W56-W63 N27-N30
      3 1103 1591 (N38W23 L=351)   381    106        1 (1075/1544-1136/1617) -90.0 W21-W25 N34-N40
      4 1657 1558 (N19W61 L= 30)  1985    324        4 (1501/1522-1782/1591) -17.9 W47-W78 N13-N23
      5  292 1390 (N45E53 L=275)  1065    146        0 ( 219/1372- 358/1410)  -5.4 E65-E43 N43-N47
      6  513 1329 (N37E28 L=301)   281     44        1 ( 488/1302- 536/1353)  33.7 E29-E25 N34-N37
      7 1760 1055 (S11W48 L= 16)   327     84        2 (1736/1016-1778/1088)  45.0 W45-W49 S12-S08
      8 1411  849 (S12W18 L=347)  1116    272        3 (1362/ 736-1440/ 987)  55.3 W13-W22 S18-S04
      9  753  720 (S04E23 L=306)   363    120        1 ( 731/ 664- 776/ 779)  56.3 E25-E20 S07-S01
     10 1090  403 (S32E13 L=316)  1378    283        3 ( 977/ 306-1154/ 503)  25.3 E24-E05 S36-S26
     11  584  429 (S19E45 L=284)   417     59        0 ( 557/ 406- 609/ 453)  18.4 E47-E41 S20-S17
     12  793  398 (S25E32 L=297)   389     99        0 ( 758/ 362- 825/ 437)  29.7 E35-E28 S27-S22
  

    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