Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1362 1643 (N43W33 L=215)   614    138        1 (1298/1617-1422/1673) -24.0 W27-W37 N41-N46
      2 1646 1498 (N31W53 L=235)   516     92        1 (1606/1473-1694/1519) -14.0 W49-W57 N28-N32
      3 1200 1484 (N33W15 L=197)   389     68        0 (1166/1460-1235/1510) -36.9 W12-W17 N31-N34
      4  231 1312 (N22E64 L=118)   600    133        0 ( 201/1269- 278/1335)  33.7 E68-E58 N19-N23
      5  506 1299 (N22E36 L=146)   551    148        0 ( 451/1257- 552/1323)  23.2 E39-E31 N19-N23
      6 1215 1249 (N17W13 L=195)   285    121        0 (1179/1197-1247/1296) -71.6 W11-W14 N14-N20
      7  360 1212 (N17E47 L=135)   348     63        2 ( 338/1187- 383/1242) -36.9 E48-E44 N15-N18
      8  844 1056 (N07E11 L=171)   210     93        0 ( 803/1044- 896/1065)   0.0 E13-E07 N05-N07
      9 1062  863 (S06W02 L=184)   692    195        1 ( 975/ 806-1131/ 923) -41.2 E03-W05 S09-S01
     10 1530  772 (S14W33 L=215)   500    113        0 (1502/ 723-1549/ 831)  81.9 W30-W34 S16-S09
     11  237  754 (S12E61 L=121)   933    157        2 ( 222/ 675- 268/ 816)  45.0 E66-E55 S17-S08
     12  712  689 (S16E21 L=161)  1295    274        4 ( 589/ 635- 829/ 754)  -7.1 E29-E13 S19-S11
     13 1735  718 (S19W52 L=234)   510    116        1 (1699/ 681-1767/ 751)  59.0 W48-W55 S21-S16
     14 1118  610 (S22W05 L=187)  1350    343        2 (1011/ 499-1250/ 715)  43.0 E02-W13 S29-S15
     15 1528  609 (S25W34 L=216)   449    118        0 (1489/ 563-1570/ 653)  50.2 W32-W37 S27-S21
     16  835  607 (S22E14 L=168)   252     53        0 ( 818/ 585- 854/ 636) -90.0 E15-E12 S23-S19
     17  983  229 (S54E09 L=173)   487     47        2 ( 962/ 200-1003/ 256) -74.1 E10-E06 S57-S50
     18 1099  206 (S57W04 L=186)   861    115        1 (1058/ 189-1157/ 234) -24.4 W00-W11 S59-S53
     19 1203  188 (S61W17 L=199)   713     69        0 (1181/ 161-1223/ 223) -83.7 W14-W21 S64-S55
  

    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