Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  827 1649 (N39E25 L=181)   356    100        0 ( 784/1621- 876/1672)  -6.3 E29-E20 N37-N40
      2 1578 1495 (N38W39 L=245)  1449    301        5 (1539/1357-1607/1622)  81.5 W34-W44 N28-N48
      3 1039 1560 (N36E06 L=200)   263     45        1 (1021/1535-1059/1585) -53.1 E07-E04 N33-N37
      4  583 1481 (N24E38 L=169)   358    159        0 ( 550/1434- 635/1507)  45.0 E39-E33 N21-N26
      5 1241 1209 (N15W11 L=218)   228    130        1 (1208/1169-1292/1263) -45.0 W09-W15 N12-N17
      6 1332  745 (S13W23 L=229)   601    277        0 (1223/ 684-1430/ 839)  47.5 W16-W28 S17-S05
      7 1178  589 (S24W16 L=222)   381    110        2 (1127/ 554-1216/ 615)  -9.5 W12-W18 S26-S22
      8  714  537 (S34E16 L=190)   243     59        0 ( 678/ 525- 750/ 546)  18.4 E19-E13 S35-S33
      9  378  495 (S43E52 L=154)  1056    129        2 ( 317/ 475- 422/ 527)   6.0 E64-E45 S45-S40
     10 1085  371 (S42W15 L=221)   304     75        0 (1044/ 358-1120/ 385)   0.0 W11-W17 S43-S40
     11 1043  232 (S56W18 L=224)   358     48        1 (1019/ 210-1063/ 254)  56.3 W15-W20 S58-S52
  

    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