Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  854 1503 (N40W00 L=127)   612    180        1 ( 832/1434- 891/1579)  73.3 E02-W03 N34-N44
      2 1082 1543 (N36W17 L=145)   311     76        1 (1057/1511-1117/1569)  11.3 W15-W20 N34-N38
      3  871 1232 (N22E05 L=122)  1033    314        0 ( 801/1117- 960/1359) -74.5 E07-E01 N13-N31
      4 1827 1137 (S06W57 L=185)  1057    229        0 (1788/1049-1862/1222) -90.0 W53-W60 S12-S01
      5  123 1141 (N28E75 L= 52)  1259     92        2 ( 102/1100- 139/1190) -59.0 E81-E72 N25-N30
      6  198  886 (N12E63 L= 64)  1931    239        6 ( 138/ 801- 236/ 944) -69.4 E72-E58 N07-N15
      7  977  858 (S01E06 L=121)   282    140        0 ( 910/ 834-1034/ 891)  -7.1 E10-E02 S02-N00
      8 1671  757 (S24W37 L=164)  1502    301        2 (1588/ 629-1733/ 876)  40.6 W29-W43 S30-S18
      9  532  740 (S00E36 L= 91)   368    103        0 ( 517/ 688- 550/ 786)  68.2 E37-E34 S03-N02
     10  967  649 (S14E11 L=116)   305    126        1 ( 931/ 603- 995/ 713) -90.0 E12-E10 S16-S08
     11 1207  676 (S17W03 L=131)   502     86        3 (1174/ 633-1246/ 707) -90.0 W01-W06 S19-S14
     12  926  548 (S19E16 L=111)   650    232        4 ( 866/ 497- 977/ 629) -84.3 E19-E13 S23-S13
     13 1085  515 (S25E08 L=120)   271    113        0 (1045/ 477-1113/ 563)  33.7 E11-E04 S26-S21
     14  398  428 (S18E62 L= 66)   815    121        0 ( 336/ 405- 439/ 460)   6.7 E72-E54 S19-S16
     15  939  295 (S38E26 L=101)  4363    662        8 ( 699/ 174-1174/ 422)  10.3 E59-E04 S43-S33
     16 1604  366 (S54W36 L=163)   845     87        2 (1571/ 341-1650/ 396)  -4.8 W30-W44 S56-S49
  

    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