Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  258 1201 (S15E56 L=177)   676    142        1 ( 216/1146- 308/1263)  58.0 E58-E52 S18-S10
      2 1790  982 (N15W50 L=282)  3195    540        8 (1730/ 766-1823/1242) -55.8 W39-W58 N03-N28
      3 1573 1120 (N15W29 L=262)   291    104        0 (1560/1061-1582/1177) -60.3 W27-W31 N11-N18
      4 1434 1072 (N08W22 L=254)   671    200        0 (1365/ 995-1515/1168) -21.0 W15-W28 N05-N11
      5  890 1042 (S09E08 L=225)   234    113        0 ( 878/ 982- 900/1099) -56.3 E09-E05 S12-S05
      6  379 1037 (S22E41 L=192)   269    121        0 ( 341/1008- 417/1069) -15.9 E44-E37 S22-S20
      7  819  980 (S14E10 L=222)   496    143        2 ( 747/ 947- 883/1007)  12.5 E14-E05 S16-S12
      8  565  986 (S20E26 L=206)   221     92        0 ( 517/ 975- 606/1000)  18.4 E29-E23 S21-S18
      9  256  773 (S41E49 L=183)   911    288        1 ( 184/ 740- 326/ 827)  31.0 E57-E42 S44-S35
     10  734  724 (S31E09 L=224)   295     98        1 ( 682/ 695- 779/ 764)  -7.1 E13-E05 S32-S29
     11 1198  673 (S20W21 L=253)   272     82        1 (1159/ 655-1248/ 699)   9.5 W17-W23 S21-S19
     12  632  499 (S49E11 L=222)  1100    253        2 ( 504/ 473- 719/ 524)  13.4 E23-E02 S52-S45
  

    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