Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1490 1544 (N33W43 L=355)   992    218        1 (1396/1500-1588/1587) -31.6 W35-W49 N29-N37
      2 1597 1371 (N21W45 L=357)   758    158        1 (1561/1292-1642/1427)  32.5 W39-W50 N16-N24
      3 1788 1266 (N11W60 L= 12)   618     65        2 (1761/1238-1816/1298) -63.4 W57-W63 N09-N13
      4 1579 1215 (N11W39 L=352)   412     85        1 (1551/1181-1608/1247) -53.1 W36-W40 N09-N13
      5 1671  652 (S25W45 L=358)   313     96        0 (1620/ 622-1706/ 675)  14.0 W40-W48 S26-S24
      6  795  580 (S21E19 L=294)   250     96        0 ( 752/ 556- 839/ 604)  15.9 E22-E15 S22-S19
      7 1086  495 (S30E00 L=312)   461    165        0 (1064/ 425-1110/ 559)  78.7 E02-W01 S34-S24
      8  862  413 (S34E18 L=295)   528    171        0 ( 840/ 337- 895/ 488)  56.3 E22-E13 S40-S28
      9  778  294 (S44E30 L=283)   653     88        1 ( 738/ 263- 818/ 326) -45.0 E33-E27 S46-S40
  

    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