Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1581 1507 (N21W50 L=330)   740    149        1 (1527/1472-1624/1525) -39.8 W45-W53 N18-N23
      2 1126 1498 (N33W19 L=298)   335     56        0 (1095/1478-1160/1523) -45.0 W16-W20 N30-N35
      3 1540 1427 (N18W43 L=323)   710    162        3 (1487/1374-1601/1464) -45.0 W40-W46 N13-N20
      4 1785 1386 (N08W64 L=344)   951    158        3 (1743/1324-1836/1453) -63.4 W61-W67 N03-N13
      5  602 1308 (N32E22 L=257)  1353    450        6 ( 534/1179- 662/1452) -70.7 E27-E18 N23-N43
      6 1744 1345 (N08W57 L=337)   746    154        2 (1702/1296-1777/1405) -74.1 W54-W59 N04-N11
      7 1283 1324 (N19W23 L=302)   490    118        2 (1237/1263-1324/1368) -82.9 W20-W24 N14-N22
      8 1748 1061 (S08W48 L=328)   891    298        1 (1680/1001-1829/1131) -42.0 W43-W56 S12-S03
      9  930 1103 (N13E04 L=276)   197     89        0 ( 895/1082- 983/1119)  -9.5 E06-W00 N12-N14
     10  183 1025 (N21E62 L=218)  2603    381        3 ( 106/ 950- 261/1088) -32.7 E70-E52 N15-N24
     11 1794  878 (S20W49 L=329)   370     89        0 (1763/ 849-1820/ 910)   8.1 W45-W52 S21-S18
     12 1382  733 (S18W16 L=295)  1315    336        4 (1296/ 615-1443/ 907)  49.1 W08-W22 S23-S08
     13  454  830 (N06E40 L=240)   325     95        1 ( 428/ 800- 476/ 864) -68.2 E41-E38 N03-N08
     14  462  640 (S06E45 L=235)   810    254        4 ( 382/ 595- 541/ 671)   4.1 E52-E38 S07-S04
     15  319  477 (S15E67 L=213)  1077    153        2 ( 290/ 440- 354/ 505) -38.7 E70-E62 S17-S13
     16  810  457 (S23E27 L=253)   492     82        2 ( 776/ 424- 841/ 488)  31.0 E29-E24 S25-S20
     17 1230  387 (S38E01 L=278)   360     41        1 (1205/ 365-1261/ 408) -45.0 E03-W00 S40-S36
  

    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