Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1204 1471 (N31W08 L=287)   315    103        0 (1190/1423-1220/1525) -81.9 W06-W09 N28-N35
      2 1380 1270 (N20W21 L=300)   254     69        1 (1353/1254-1421/1287) -11.3 W18-W23 N18-N20
      3 1781 1118 (N14W54 L=333)   716    177        0 (1740/1048-1808/1190) -60.9 W50-W56 N09-N18
      4  191 1003 (S09E63 L=217)   548    107        0 ( 145/ 973- 240/1022)  -8.1 E70-E56 S10-S07
      5  142  925 (S14E71 L=208)  1069     92        1 ( 106/ 886- 178/ 952) -13.0 E82-E65 S16-S12
      6  440  867 (S14E38 L=242)   387     72        1 ( 404/ 842- 462/ 904)  26.6 E40-E35 S15-S11
      7  999  855 (S09W00 L=280)   285     87        1 ( 969/ 817-1031/ 891)  90.0 E01-W01 S11-S06
      8  296  779 (S22E52 L=228)  1443    217        3 ( 253/ 709- 359/ 870) -45.0 E55-E45 S26-S16
      9  498  814 (S17E33 L=246)   299     65        0 ( 486/ 770- 513/ 851)  78.7 E33-E31 S19-S14
     10 1691  578 (S20W57 L=336)   463     73        1 (1670/ 539-1715/ 610)  78.7 W54-W58 S23-S18
     11 1612  420 (S32W59 L=338)   549     95        1 (1567/ 388-1638/ 441)  -9.5 W52-W64 S34-S30
     12 1383  349 (S40W40 L=319)   445     47        0 (1366/ 320-1399/ 377) -51.3 W37-W43 S42-S37
  

    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