Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  935 1609 (N37E16 L=356)   501    200        0 ( 877/1553-1012/1673) -26.6 E23-E09 N33-N40
      2 1733 1445 (N37W54 L= 67)   533     49        0 (1701/1428-1763/1465)   7.1 W50-W59 N36-N38
      3  494 1385 (N15E41 L=331)   599    189        1 ( 459/1296- 516/1460) -45.0 E46-E37 N10-N19
      4  229 1371 (N11E68 L=304)   777     69        2 ( 209/1336- 245/1414) -59.0 E71-E64 N08-N13
      5 1647 1347 (N29W41 L= 54)   390    108        1 (1608/1307-1685/1385) -39.8 W38-W44 N26-N31
      6  738 1233 (N09E21 L=352)   851    205        0 ( 693/1137- 779/1323)  80.5 E22-E18 N03-N15
      7  952 1196 (N10E07 L=  5)   231    103        1 ( 920/1148- 977/1238) -51.3 E09-E04 N07-N12
      8  423 1119 (S02E41 L=331)   379     87        0 ( 396/1092- 456/1156) -21.8 E43-E37 S03-N00
      9  950  665 (S23W00 L= 12)  1903    652        4 ( 699/ 569-1224/ 771)  32.8 E16-W16 S32-S12
     10  323  607 (S36E53 L=319)   925    202        0 ( 279/ 530- 388/ 696) -57.5 E56-E47 S41-S30
     11  462  635 (S32E36 L=336)   444    114        0 ( 401/ 621- 517/ 651)   5.2 E42-E31 S32-S31
  

    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