Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  736 1518 (N40E16 L=336)   231     51        0 ( 724/1487- 752/1543)  63.4 E17-E14 N38-N42
      2 1415 1501 (N30W35 L= 26)  3913    887       11 (1122/1441-1755/1541) -13.2 W13-W64 N22-N36
      3 1627 1308 (N15W45 L= 37)   303    107        0 (1582/1293-1667/1323) -23.2 W41-W48 N13-N16
      4  183 1305 (N29E65 L=287)   802    238        1 ( 147/1280- 217/1322)  -6.3 E69-E60 N27-N29
      5  904 1267 (N22E05 L=347)   239     65        1 ( 877/1232- 929/1292) -63.4 E06-E03 N19-N23
      6 1020 1257 (N20W02 L=354)   192     87        0 ( 969/1237-1055/1286)  21.8 E01-W04 N19-N21
      7 1545  340 (S51W40 L= 31)  1040    163        2 (1511/ 307-1577/ 383) -45.0 W34-W46 S55-S45
      8 1004  189 (S56E18 L=334)   897    125        0 ( 945/ 171-1061/ 212) -37.9 E22-E12 S59-S52
  

    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