Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1020 1510 (N36W02 L=121)   964    286        1 ( 918/1426-1183/1591) -33.7 E05-W13 N29-N42
      2 1450 1329 (N21W30 L=149)   484    134        0 (1401/1285-1509/1370)  24.0 W26-W35 N18-N23
      3  504 1336 (N25E36 L= 83)   743    137        2 ( 441/1305- 559/1368)  18.4 E40-E31 N23-N27
      4 1316 1308 (N21W21 L=139)   297     70        3 (1292/1285-1344/1331) -45.0 W18-W22 N18-N22
      5 1685 1270 (N15W48 L=167)   920    170        4 (1638/1198-1726/1319)  33.7 W43-W52 N11-N18
      6 1679 1110 (N06W45 L=164)   517     56        1 (1654/1073-1704/1143) -53.1 W43-W47 N03-N07
      7 1438  614 (S24W27 L=145)  1115    327        4 (1386/ 483-1488/ 731)  76.8 W24-W29 S33-S16
      8  740  410 (S36E25 L= 94)   886    252        1 ( 666/ 354- 819/ 477)  29.4 E33-E17 S40-S30
  

    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