Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  837 1650 (N42E21 L=286)  3933    650       11 ( 558/1512-1115/1760) -11.7 E56-W02 N33-N47
      2 1345 1530 (N37W21 L=328)  1388    253        3 (1226/1503-1466/1555)   3.0 W11-W30 N34-N38
      3 1580 1434 (N31W39 L=347)   932    170        2 (1521/1364-1641/1485) -48.8 W34-W43 N26-N34
      4  325 1376 (N18E55 L=253)  1140    117        4 ( 268/1339- 376/1420)  21.8 E60-E50 N15-N20
      5  572 1369 (N20E33 L=274)   590     93        2 ( 547/1326- 609/1419) -63.4 E35-E30 N17-N23
      6 1405 1316 (N22W23 L=331)   326    106        0 (1391/1268-1416/1376)  81.9 W22-W24 N19-N26
      7 1318 1309 (N21W17 L=325)   223     76        0 (1295/1277-1345/1346) -45.0 W15-W19 N19-N23
      8 1807 1244 (N20W58 L=  6)   715    147        0 (1734/1223-1847/1282)   8.1 W50-W64 N18-N22
      9  176 1222 (N07E68 L=239)   655     67        0 ( 153/1204- 213/1250)  -5.7 E72-E62 N06-N08
     10  526 1108 (N03E32 L=275)  1233    250        1 ( 497/ 994- 568/1220)  90.0 E34-E29 S03-N10
     11 1844 1060 (N09W61 L=  8)   456     51        0 (1821/1033-1867/1086) -26.6 W57-W63 N06-N10
     12  390  997 (S04E42 L=266)   335     65        1 ( 357/ 977- 431/1013)  18.4 E44-E38 S05-S03
     13 1832  926 (N00W60 L=  7)   565     96        0 (1797/ 904-1873/ 963)  14.0 W56-W64 S01-N02
     14  220  799 (S18E60 L=247)   469     55        1 ( 206/ 769- 238/ 833)  53.1 E62-E57 S20-S16
     15  759  737 (S18E15 L=292)   512    183        1 ( 673/ 706- 836/ 782) -15.3 E20-E09 S19-S15
     16 1412  577 (S24W30 L=338)  3118    593        6 (1161/ 472-1656/ 688)  25.9 W14-W47 S32-S15
     17  542  659 (S24E31 L=277)   356     86        2 ( 506/ 637- 584/ 687)  29.7 E34-E27 S26-S22
     18  358  535 (S35E53 L=255)   387     53        0 ( 324/ 523- 387/ 546)   7.1 E57-E48 S36-S34
     19 1040  436 (S36W06 L=313)   462    101        2 ( 993/ 407-1078/ 479)  56.3 W02-W08 S38-S32
  

    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