Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  840 1847 (N65W26 L=293)   613     65        1 ( 817/1813- 859/1879)  19.7 W20-W34 N61-N67
      2 1359 1854 (N42W67 L=334)   794     78        1 (1316/1827-1402/1876) -11.3 W60-W76 N40-N43
      3  514 1661 (N61E25 L=242)  1267    244        2 ( 417/1603- 614/1723)   5.9 E39-E09 N58-N63
      4 1438 1437 (N17W37 L=303)   888    268        1 (1312/1381-1553/1481) -45.0 W30-W42 N10-N22
      5  869 1331 (N27E01 L=266)   828    256        3 ( 760/1291- 991/1367) -42.3 E06-W06 N21-N31
      6 1811 1198 (S07W56 L=322)   480     90        1 (1794/1153-1843/1231)  26.6 W52-W60 S09-S05
      7 1705 1073 (S11W42 L=309)   620    135        0 (1679/1007-1723/1139)  50.2 W38-W45 S13-S06
      8  685  987 (N12E20 L=246)   201    111        0 ( 640/ 961- 724/1004) -45.0 E23-E18 N09-N13
      9 1717  742 (S30W38 L=304)   322     57        1 (1696/ 714-1738/ 775)  45.0 W35-W39 S32-S28
     10  185  723 (N07E69 L=197)   478     35        0 ( 160/ 710- 212/ 736)   0.0 E73-E65 N06-N08
     11 1523  565 (S35W19 L=286)  1481    253        5 (1414/ 521-1620/ 635)   0.0 W10-W27 S37-S31
     12  324  560 (S04E61 L=206)   393     73        1 ( 300/ 517- 344/ 596)  71.6 E63-E58 S07-S01
     13  675  493 (S15E36 L=231)   424    117        4 ( 634/ 463- 709/ 523)   0.0 E39-E33 S17-S13
     14  452  408 (S16E58 L=209)  1552    243        2 ( 402/ 336- 518/ 498)  79.7 E62-E50 S21-S10
     15  560  369 (S21E51 L=215)  1052    142        3 ( 514/ 316- 616/ 428) -83.7 E54-E48 S25-S16
     16  837  357 (S28E32 L=234)   316    100        1 ( 793/ 343- 871/ 383) -45.0 E34-E28 S29-S25
     17  686  316 (S27E45 L=221)   553     83        2 ( 650/ 279- 715/ 360)  39.8 E48-E42 S29-S24
     18 1003  229 (S43E31 L=236)  2767    457        4 ( 812/ 129-1149/ 333)   5.7 E55-E14 S47-S38
  

    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