Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1491 1448 (N28W37 L= 15)   306     92        0 (1465/1408-1514/1491) -90.0 W36-W38 N24-N30
      2 1730 1362 (N20W56 L= 34)   554     93        1 (1700/1331-1768/1407) -51.3 W53-W59 N17-N22
      3  219 1250 (N20E63 L=276)   719    105        2 ( 196/1208- 255/1297) -33.7 E67-E57 N17-N23
      4  582  985 (N04E28 L=310)   824    231        1 ( 547/ 887- 611/1095)  79.7 E31-E26 S01-N10
      5  279  891 (S01E53 L=285)  2987    492        1 ( 224/ 727- 421/1067) -55.0 E58-E42 S11-N09
      6 1033  866 (S05E00 L=338)   329     93        0 ( 991/ 844-1075/ 899)  26.6 E02-W02 S06-S03
      7 1252  734 (S15W13 L=351)   443    122        2 (1190/ 699-1293/ 762) -33.7 W08-W15 S17-S12
      8  475  687 (S14E39 L=299)   312     57        1 ( 438/ 667- 502/ 702)   0.0 E42-E36 S15-S13
      9  469  580 (S21E42 L=296)   788    186        0 ( 386/ 540- 549/ 625) -22.6 E48-E36 S24-S18
     10  662  343 (S40E36 L=303)  1658    281        1 ( 572/ 247- 752/ 449)  25.9 E57-E23 S49-S31
  

    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