Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1252 1813 (N51E12 L=276)   528     57        1 (1224/1791-1280/1839)  36.9 E14-E09 N48-N54
      2  472 1709 (N21E60 L=228)   631    116        2 ( 444/1683- 503/1726)  11.3 E63-E56 N19-N22
      3  665 1540 (N14E36 L=252)   659    294        0 ( 636/1435- 701/1653) -42.7 E42-E29 N08-N20
      4  301 1202 (S13E51 L=237)   589     96        1 ( 279/1148- 340/1245) -63.4 E54-E48 S16-S10
      5  121 1169 (S17E71 L=217)  1421    171        1 ( 100/1110- 140/1244) -53.1 E73-E66 S20-S12
      6  181 1186 (S16E63 L=225)  1527    284        8 ( 147/1111- 221/1235) -82.9 E68-E58 S20-S12
      7  226 1120 (S19E56 L=231)   714    157        4 ( 184/1085- 268/1151)  51.3 E60-E52 S21-S16
      8  746 1089 (S10E17 L=271)   249     53        1 ( 714/1066- 777/1108)  45.0 E19-E15 S11-S08
      9  514 1067 (S17E32 L=256)   590    233        1 ( 427/1039- 639/1094)  15.9 E37-E23 S19-S14
     10 1163  976 (S06W09 L=297)   402    120        2 (1113/ 922-1210/1024) -15.9 W05-W12 S08-S04
     11  181  979 (S28E58 L=230)  1287    483        1 (  84/ 941- 278/1005)  14.0 E68-E48 S30-S25
     12  649  886 (S24E19 L=269)  2343    749        1 ( 345/ 802- 977/ 977)   6.5 E42-W03 S31-S16
     13 1188  690 (S21W19 L=307)   290     76        1 (1149/ 657-1229/ 720)  -8.1 W15-W22 S21-S20
     14  875  583 (S36W04 L=292)   742    198        1 ( 776/ 541- 973/ 612)   8.1 E03-W11 S38-S34
     15 1344  560 (S23W33 L=321)   269     83        1 (1324/ 526-1370/ 599) -31.0 W31-W36 S24-S21
     16 1721  531 (S12W63 L=351)  2495    404       10 (1654/ 449-1803/ 595)  47.7 W57-W71 S17-S06
  

    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