Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  763 1760 (N30E42 L= 61)   330     46        0 ( 737/1740- 791/1782)  -9.5 E45-E39 N29-N31
      2  495 1673 (N19E57 L= 46)   554     78        1 ( 477/1635- 513/1711) -50.2 E62-E53 N15-N21
      3 1107 1386 (N15E05 L= 98)   562    202        0 (1065/1292-1150/1473) -45.0 E09-W00 N10-N19
      4 1234 1349 (N17W03 L=106)   435     96        1 (1189/1329-1279/1372)  11.3 W00-W05 N15-N19
      5 1336 1254 (N14W11 L=114)   253     97        0 (1319/1205-1361/1301) -45.0 W09-W14 N11-N16
      6 1746  774 (N01W53 L=156)   649     68        1 (1724/ 737-1766/ 818) -18.4 W50-W56 S01-N02
      7 1718  651 (S06W56 L=159)   339     39        1 (1699/ 628-1739/ 670)   0.0 W54-W58 S07-S05
      8 1068  532 (S33W18 L=121)  2064    513        2 ( 890/ 473-1283/ 655)  49.4 W08-W26 S42-S20
      9 1459  500 (S23W45 L=148)   443     57        1 (1430/ 474-1489/ 526)   0.0 W41-W48 S24-S21
     10 1546  338 (S28W66 L=169)   763     78        3 (1522/ 300-1572/ 374) -45.0 W61-W69 S30-S25
     11 1117  297 (S46W37 L=140)   389     95        0 (1066/ 286-1170/ 309)  39.8 W32-W40 S48-S43
  

    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