Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  717 1845 (N36E53 L=169)   536    118        0 ( 698/1816- 736/1874) -26.6 E57-E49 N33-N38
      2  929 1268 (N05E12 L=210)   396    127        1 ( 855/1245- 992/1290)  14.0 E16-E08 N02-N05
      3  358 1156 (S15E45 L=176)   542    128        0 ( 315/1102- 404/1197)  60.3 E48-E42 S19-S12
      4 1846 1156 (N27W51 L=273)   503     70        1 (1823/1126-1867/1188) -56.3 W48-W53 N25-N29
      5 1690 1093 (N18W38 L=260)   264     63        0 (1669/1055-1705/1126) -45.0 W36-W40 N15-N19
      6 1569  991 (N08W32 L=254)   860    183        0 (1511/ 912-1628/1068) -24.4 W26-W38 N04-N10
      7  477  998 (S22E33 L=189)   320    104        1 ( 447/ 962- 519/1050) -26.6 E35-E29 S23-S19
      8 1883  896 (N14W63 L=284)  3065    353        2 (1847/ 727-1902/1027) -50.9 W54-W72 N05-N21
      9  694  923 (S21E16 L=206)   272     96        0 ( 646/ 910- 746/ 938)  15.9 E19-E12 S21-S19
     10  339  744 (S41E41 L=181)  1044    317        1 ( 243/ 717- 455/ 769)  16.7 E50-E30 S44-S36
     11 1016  614 (S29W12 L=234)   330    119        0 ( 953/ 602-1079/ 629)  29.7 W08-W15 S30-S26
     12  793  437 (S48W05 L=227)   335     41        1 ( 760/ 413- 820/ 462)   0.0 W01-W07 S49-S46
  

    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