Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1428 1701 (N35W53 L=296)  1152    206        2 (1344/1660-1521/1722) -36.0 W46-W61 N30-N38
      2 1628 1696 (N28W74 L=317)   673     81        1 (1595/1679-1669/1713) -14.0 W68-W80 N25-N29
      3  404 1477 (N49E42 L=201)   475     79        1 ( 379/1449- 437/1512) -26.6 E46-E36 N46-N51
      4 1662 1386 (N11W51 L=294)  1074    448        4 (1486/1273-1777/1465) -51.7 W41-W57 N01-N20
      5 1389 1275 (N13W28 L=270)   905    276        1 (1348/1144-1446/1381)  86.2 W25-W31 N04-N19
      6  755 1072 (N16E15 L=227)   917    480        6 ( 639/ 891- 871/1222) -71.6 E20-E10 N02-N26
      7 1024  891 (S00E03 L=239)  1036    416        2 ( 918/ 746-1171/ 942)  10.0 E12-W05 S06-N03
      8  134  908 (N15E73 L=169)   831     65        2 ( 110/ 882- 158/ 939) -56.3 E78-E69 N13-N16
      9 1814  814 (S26W51 L=293)   744     76        1 (1782/ 776-1836/ 856) -80.5 W47-W54 S29-S23
     10  808  723 (S06E19 L=223)   325    119        2 ( 763/ 679- 842/ 765) -90.0 E21-E17 S08-S03
     11  280  553 (S08E67 L=175)  5514    529        8 ( 137/ 393- 386/ 734) -67.8 E82-E62 S19-N04
     12  372  553 (S09E56 L=186)  2284    445        6 ( 262/ 394- 445/ 717)  84.8 E62-E51 S20-N02
     13 1736  595 (S39W44 L=286)  1150    188        3 (1686/ 538-1795/ 686)  32.5 W38-W50 S42-S34
     14  508  589 (S08E43 L=199)   731    140        1 ( 482/ 543- 558/ 651) -81.9 E45-E40 S11-S04
     15  885  503 (S21E21 L=221)   625    126        2 ( 853/ 445- 910/ 564)  54.5 E23-E17 S24-S17
     16  754  480 (S20E30 L=212)   365    108        3 ( 716/ 448- 792/ 518)  23.2 E33-E26 S21-S17
     17 1630  480 (S44W33 L=275)   704    138        1 (1565/ 452-1680/ 511)  -8.7 W26-W39 S46-S42
     18  523  362 (S24E55 L=188)   576    117        1 ( 482/ 308- 556/ 400)  21.8 E64-E49 S27-S20
  

    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