Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  386 1653 (N19E64 L= 30)   725     77        2 ( 358/1617- 417/1689) -63.4 E66-E59 N16-N22
      2 1682 1466 (N36W36 L=130)   454     85        1 (1648/1437-1722/1483)  32.0 W33-W41 N33-N38
      3 1421 1447 (N27W15 L=108)   246     83        0 (1386/1410-1454/1480)  63.4 W13-W17 N23-N29
      4 1765 1434 (N37W46 L=139)   446     50        0 (1737/1411-1788/1456)  -9.5 W42-W48 N35-N38
      5  155 1128 (S15E66 L= 27)  4054    389       10 ( 117/ 936- 216/1239)  66.0 E74-E60 S26-S08
      6 1266 1178 (N07W10 L=104)   352    136        1 (1207/1161-1335/1200)   7.1 W06-W14 N05-N08
      7  479 1103 (S13E35 L= 58)   361    127        1 ( 420/1064- 526/1125) -14.0 E39-E30 S14-S11
      8 1084  897 (S13W06 L= 99)  1024    532        0 ( 865/ 821-1290/ 977)  37.3 E05-W17 S21-S05
      9  486  925 (S23E32 L= 61)   434    126        2 ( 417/ 904- 544/ 953)  21.8 E37-E27 S24-S20
     10 1209  642 (S24W20 L=114)   337    102        0 (1166/ 628-1264/ 662)  18.4 W17-W23 S26-S22
     11 1350  523 (S28W34 L=127)   284     49        1 (1325/ 496-1371/ 545) -14.0 W32-W36 S28-S25
  

    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