Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1315 1753 (N39W48 L=287)  1891    246        4 (1204/1716-1445/1778) -19.3 W36-W58 N33-N42
      2 1607 1576 (N19W57 L=296)  4322    717        6 (1458/1437-1754/1682) -67.2 W49-W69 N09-N28
      3  776 1377 (N30E06 L=233)   553    198        3 ( 691/1344- 869/1410) -18.4 E12-W00 N26-N31
      4  522 1369 (N36E25 L=214)   430    143        1 ( 464/1331- 594/1397)   0.0 E30-E19 N34-N36
      5 1108 1178 (N10W08 L=248)   457    153        1 (1063/1109-1149/1241) -83.7 W07-W09 N05-N14
      6 1813 1176 (S08W56 L=295)   590    105        0 (1784/1122-1839/1219)  33.7 W51-W60 S10-S05
      7 1643 1040 (S11W38 L=277)   312     73        0 (1620/1006-1666/1079) -90.0 W36-W39 S12-S07
      8 1862  880 (S27W57 L=297)   700    166        0 (1825/ 802-1894/ 954)  41.2 W52-W62 S30-S22
      9  641  775 (S01E28 L=211)   211     49        1 ( 620/ 746- 657/ 798)  56.3 E29-E26 S03-N00
     10 1738  721 (S33W42 L=282)   713    125        1 (1706/ 680-1778/ 783) -48.4 W39-W47 S36-S27
     11 1586  674 (S31W28 L=267)   857    241        1 (1514/ 603-1661/ 776)  26.6 W21-W35 S34-S26
     12  188  726 (N04E69 L=170)   614     78        0 ( 172/ 692- 205/ 774) -90.0 E71-E66 N02-N07
     13  743  706 (S08E24 L=215)   444    157        0 ( 701/ 627- 763/ 768)  83.7 E25-E22 S12-S03
     14  638  655 (S08E32 L=207)   220     32        0 ( 613/ 641- 663/ 672) -33.7 E33-E30 S09-S07
     15 1014  566 (S22E12 L=227)   625    149        1 ( 977/ 515-1034/ 640)  36.9 E16-E08 S25-S18
     16 1440  456 (S42W13 L=252)  1226    230        1 (1336/ 394-1509/ 526)  13.2 W03-W20 S44-S38
     17 1045  204 (S50E29 L=210)  1866    272        4 ( 911/ 161-1156/ 261)  -2.1 E43-E16 S52-S47
  

    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