Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  839 1702 (N51W10 L=279)   315     60        0 ( 804/1689- 872/1715) -31.0 W07-W13 N49-N52
      2 1128 1409 (N23W17 L=286)   509    128        0 (1098/1358-1185/1452)  20.6 W13-W21 N20-N24
      3  984 1369 (N25W07 L=276)   408     91        1 ( 942/1336-1027/1396) -21.8 W04-W09 N23-N25
      4  208 1334 (N42E61 L=208)   616     74        1 ( 165/1322- 249/1347)  -6.7 E71-E54 N41-N43
      5 1715  888 (S22W41 L=310)   708    178        0 (1680/ 817-1737/ 972)  60.9 W36-W42 S25-S16
      6 1509  875 (S17W25 L=294)   581    187        0 (1450/ 845-1567/ 927) -58.0 W21-W29 S19-S11
      7  318  778 (N06E52 L=217)   637     86        2 ( 282/ 749- 363/ 808)   0.0 E56-E47 N04-N08
      8  326  659 (S00E56 L=214)   662     72        1 ( 298/ 620- 355/ 695)  31.0 E59-E51 S02-N01
      9  977  523 (S23E16 L=253)   766    280        1 ( 857/ 486-1069/ 592) -12.1 E23-E09 S25-S19
     10 1291  473 (S36W01 L=270)   336     53        1 (1264/ 450-1316/ 493)  18.4 W00-W03 S37-S33
     11  724  384 (S26E38 L=231)   478    125        0 ( 706/ 321- 744/ 441)  53.1 E41-E34 S30-S22
  

    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