Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1677 1543 (N34W56 L= 86)   767     85        1 (1640/1506-1708/1571) -53.1 W52-W57 N31-N36
      2  530 1480 (N34E37 L=354)   927    203        1 ( 457/1430- 614/1519) -14.9 E44-E29 N30-N36
      3  366 1285 (N20E47 L=344)  2806    700        7 ( 230/1083- 465/1442) -38.9 E64-E36 N08-N29
      4  934 1385 (N27E05 L= 25)   353     99        1 ( 870/1367- 984/1399)   7.1 E09-E01 N25-N27
      5 1792 1225 (N13W58 L= 88)   467    123        1 (1763/1174-1810/1298) -82.9 W55-W59 N09-N17
      6 1151 1181 (N13W08 L= 39)   295    106        0 (1095/1143-1204/1228) -35.5 W04-W11 N10-N16
      7 1549  987 (S00W34 L= 65)   295     83        1 (1508/ 970-1594/ 999)   0.0 W31-W37 N00-N00
      8  485  915 (S02E36 L=355)   303     66        0 ( 454/ 901- 520/ 930) -11.3 E37-E32 S03-S01
      9  770  833 (S07E16 L= 14)   384    106        0 ( 710/ 804- 821/ 864) -23.2 E20-E13 S09-S05
     10  355  647 (S20E50 L=341)   616    143        0 ( 332/ 599- 376/ 691)  49.4 E53-E47 S23-S16
     11 1424  579 (S26W27 L= 58)   555    121        1 (1359/ 565-1473/ 598)   0.0 W22-W30 S26-S24
     12  576  523 (S28E34 L=357)   823    181        3 ( 526/ 490- 641/ 559)  24.4 E38-E27 S30-S25
  

    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