Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1513 1723 (N31W61 L=233)  1114    188        1 (1411/1700-1585/1747) -20.6 W51-W68 N27-N34
      2 1158 1673 (N38W29 L=201)  1168    269        0 (1034/1622-1283/1715)  -5.2 W18-W40 N36-N39
      3 1347 1452 (N20W31 L=202)   252     85        0 (1303/1431-1380/1476) -51.3 W28-W32 N17-N22
      4  377 1288 (N33E39 L=133)   872    156        1 ( 296/1270- 457/1303) -13.2 E48-E31 N30-N35
      5  899 1199 (N16E03 L=168)   233     84        0 ( 857/1186- 945/1212) -31.0 E05-W00 N14-N17
      6  540 1162 (N22E27 L=145)   412    118        1 ( 490/1118- 596/1186)   0.0 E30-E22 N19-N22
      7 1785 1039 (S14W51 L=222)  1243    199        3 (1721/ 964-1850/1130)  23.2 W43-W59 S18-S10
      8 1651 1004 (S13W38 L=210)   565    273        1 (1564/ 952-1709/1083)  26.6 W31-W44 S15-S09
      9 1552 1057 (S07W32 L=204)   256     68        2 (1515/1033-1576/1078) -11.3 W29-W34 S08-S05
     10  183  970 (N17E62 L=110)  2176    228        3 (  90/ 922- 253/1044)   6.3 E80-E52 N14-N21
     11 1848  896 (S25W56 L=228)   802     81        2 (1815/ 865-1875/ 942)  53.1 W51-W59 S26-S21
     12 1281  869 (S11W11 L=183)   194     22        0 (1264/ 851-1298/ 891)  63.4 W10-W12 S12-S10
     13 1418  740 (S22W17 L=189)   260     58        0 (1404/ 703-1432/ 776)  63.4 W15-W18 S24-S20
     14 1224  638 (S24W03 L=175)   682    244        2 (1171/ 531-1282/ 746)  47.7 E02-W08 S29-S18
     15 1029  621 (S20E09 L=162)   726    266        2 ( 983/ 509-1082/ 730)  42.5 E15-E03 S26-S13
     16  710  483 (S21E34 L=138)  1033    385        0 ( 656/ 325- 770/ 656)  38.0 E46-E23 S30-S11
     17 1278  545 (S31W04 L=176)   394    128        0 (1215/ 518-1326/ 583)   0.0 W00-W08 S32-S29
     18 1539  412 (S50W24 L=196)   857    133        0 (1474/ 391-1593/ 447)   0.0 W17-W30 S51-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