Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1652 1459 (N24W46 L=243)   333     25        1 (1631/1439-1672/1482) -71.6 W44-W48 N23-N26
      2  322 1323 (N16E50 L=147)   314     68        1 ( 304/1301- 354/1351) -31.0 E51-E46 N14-N18
      3  657 1313 (N15E23 L=173)   684    183        2 ( 582/1289- 748/1335)   0.0 E28-E17 N13-N16
      4  272 1234 (N11E53 L=143)   976    150        3 ( 252/1160- 301/1300) -83.7 E55-E50 N06-N15
      5  209 1216 (N10E60 L=137)   778    109        0 ( 184/1162- 234/1261)  56.3 E62-E57 N07-N13
      6  310 1151 (N06E48 L=148)  1063    209        2 ( 246/1054- 362/1227)  48.0 E54-E44 N00-N10
      7  875  836 (S14E09 L=187)   936    315        0 ( 810/ 699- 910/ 982)  67.6 E14-E07 S22-S05
      8 1354  814 (S16W21 L=217)   212     73        0 (1334/ 769-1376/ 848) -59.0 W19-W22 S18-S13
      9 1430  792 (S17W26 L=223)   272    112        0 (1379/ 777-1475/ 812)   0.0 W22-W29 S17-S15
     10  881  473 (S38E11 L=185)   745    219        1 ( 770/ 453- 975/ 507)  10.0 E20-E03 S39-S35
     11 1690  442 (S39W63 L=260)   693     91        1 (1647/ 423-1733/ 464)  -6.0 W55-W74 S40-S38
     12  896  317 (S51E13 L=184)   303     57        1 ( 871/ 299- 922/ 341)  11.3 E15-E10 S52-S48
  

    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