Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1373 1212 (N12W25 L=  2)   401    125        0 (1348/1156-1404/1258) -81.9 W23-W26 N08-N15
      2  257 1166 (N22E56 L=282)   336     36        1 ( 236/1144- 281/1184) -33.7 E59-E53 N21-N23
      3 1625  889 (S11W37 L= 15)   336     88        0 (1592/ 855-1650/ 926)  31.0 W34-W39 S13-S09
      4 1855  755 (S26W63 L= 41)   540     74        1 (1829/ 736-1878/ 790)  31.0 W59-W67 S27-S23
      5 1535  657 (S25W30 L=  7)   703    204        1 (1478/ 573-1588/ 745)  51.3 W25-W33 S29-S19
      6 1372  661 (S21W17 L=355)   378    115        0 (1348/ 614-1393/ 721)  66.8 W15-W19 S24-S17
      7 1041  467 (S30E07 L=330)   403     97        0 (1021/ 415-1072/ 509)  50.2 E09-E04 S33-S27
      8  380  424 (S26E65 L=273)  1254    116        2 ( 365/ 365- 411/ 478)  26.6 E78-E58 S31-S21
      9  800  293 (S40E32 L=306)   370     47        0 ( 775/ 272- 829/ 313)  15.9 E35-E28 S41-S38
  

    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