Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1652 1424 (N38W32 L=166)  1527    289        0 (1539/1329-1734/1480)  -5.7 W21-W41 N34-N41
      2 1474 1366 (N28W18 L=152)   752    118        2 (1417/1316-1516/1411) -29.7 W13-W21 N25-N29
      3 1333 1350 (N22W09 L=143)   303     98        0 (1304/1301-1352/1399)  90.0 W08-W09 N18-N25
      4  385 1342 (S02E50 L= 84)  1213    205        4 ( 343/1275- 443/1387) -26.6 E55-E44 S05-N00
      5 1217 1256 (N13W05 L=138)   297    115        0 (1159/1235-1268/1280)  39.8 W01-W07 N10-N16
      6 1741 1091 (N20W44 L=178)   460     81        2 (1715/1048-1761/1127) -59.0 W42-W47 N17-N22
      7 1513  906 (N02W32 L=165)   758    239        2 (1488/ 807-1538/1028) -71.6 W28-W34 S03-N09
      8 1864  943 (N16W61 L=195)   585     69        3 (1839/ 916-1894/ 969)  39.8 W57-W65 N13-N18
      9 1260  890 (S06W17 L=150)  1237    324        2 (1137/ 793-1335/ 944)   4.1 W10-W24 S09-S01
     10 1663  906 (N07W42 L=176)   295     71        0 (1625/ 892-1700/ 919)  18.4 W39-W45 N05-N08
     11  983  862 (S15W02 L=136)   195    100        0 ( 930/ 841-1028/ 876)  51.3 W00-W04 S17-S12
     12  904  722 (S25W01 L=135)   360    147        0 ( 879/ 667- 924/ 796) -63.4 E01-W03 S28-S20
     13  544  779 (S31E24 L=109)   253     64        0 ( 502/ 759- 571/ 795)  38.7 E27-E21 S33-S29
     14  503  685 (S39E26 L=107)   286     37        1 ( 484/ 661- 522/ 712)  76.0 E27-E24 S40-S36
     15 1702  511 (S12W65 L=199)  1416    213        2 (1635/ 452-1766/ 591)  74.7 W61-W72 S17-S06
     16 1581  519 (S16W53 L=187)   439    100        0 (1548/ 503-1631/ 539)  23.2 W49-W57 S17-S14
     17 1483  326 (S30W61 L=195)  1673    167        4 (1422/ 283-1560/ 371)   2.4 W52-W76 S31-S27
     18 1148  191 (S51W53 L=187)   430     58        1 (1137/ 158-1161/ 227) -21.0 W47-W60 S53-S48
  

    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