Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1294 1674 (N35W37 L=133)  1034    189        1 (1217/1647-1387/1697) -17.1 W30-W43 N31-N37
      2 1557 1595 (N23W53 L=148)   474     94        0 (1520/1558-1596/1620)   5.7 W47-W57 N22-N24
      3 1104 1336 (N19W12 L=107)  1391    413        1 ( 922/1275-1270/1443)   0.0 W00-W24 N16-N21
      4  577 1377 (N33E22 L= 74)  1839    328        3 ( 451/1272- 695/1441) -33.7 E32-E14 N24-N38
      5 1844 1217 (S05W62 L=157)   404     48        0 (1827/1184-1865/1235)  21.8 W59-W64 S06-S03
      6 1731  881 (S21W44 L=139)  1486    339        1 (1666/ 755-1806/1046)  43.0 W37-W53 S28-S13
      7 1882  840 (S28W62 L=158)   565     60        1 (1858/ 809-1906/ 866)  18.4 W58-W66 S29-S26
      8  861  716 (S12E16 L= 80)   529    178        1 ( 843/ 642- 888/ 821)  56.3 E18-E12 S15-S05
      9 1386  655 (S27W15 L=110)   501     93        0 (1354/ 601-1411/ 702)  45.0 W11-W16 S30-S24
     10 1325  521 (S34W08 L=104)   228     56        0 (1292/ 507-1353/ 538) -11.3 W05-W10 S35-S32
  

    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