Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1694 1427 (N22W56 L=189)   364     45        2 (1669/1403-1717/1446) -45.0 W53-W58 N20-N23
      2 1805 1381 (N17W67 L=200)   953     88        0 (1762/1359-1847/1400)  -8.7 W62-W75 N15-N19
      3 1869 1112 (N00W65 L=198)   688     68        1 (1840/1077-1888/1148)  45.0 W61-W69 S01-N02
      4 1170 1070 (N07W09 L=142)   457    101        2 (1144/1033-1212/1125) -90.0 W07-W11 N04-N10
      5 1362  801 (S12W19 L=152)   336    161        0 (1336/ 743-1392/ 853)  66.8 W17-W21 S15-S08
      6 1759  735 (S21W52 L=185)  1715    242        2 (1673/ 666-1835/ 834)  28.1 W44-W60 S24-S16
      7 1169  743 (S13W07 L=139)   380    159        0 (1094/ 721-1220/ 773)   0.0 W02-W10 S14-S11
      8 1564  509 (S35W37 L=170)   551     80        2 (1531/ 479-1602/ 544)   0.0 W33-W40 S37-S32
      9  404  516 (S24E53 L= 80)   437     73        1 ( 371/ 498- 438/ 537)   0.0 E57-E49 S24-S22
     10 1381  314 (S49W25 L=157)   370     57        1 (1359/ 296-1412/ 339) -15.9 W21-W28 S50-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