Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1608 (N37E17 L=356)   513    165        0 ( 875/1556- 999/1667) -26.6 E22-E10 N33-N40
      2 1255 1550 (N37W09 L= 21)   258     36        1 (1231/1535-1283/1568)  14.0 W07-W11 N36-N38
      3 1728 1446 (N37W54 L= 66)   480     33        0 (1700/1427-1752/1466)   8.1 W49-W57 N35-N38
      4  488 1383 (N15E42 L=331)   716    202        0 ( 456/1289- 509/1464) -48.0 E47-E37 N09-N19
      5  224 1370 (N11E69 L=304)   894     87        1 ( 201/1330- 240/1416) -51.3 E72-E65 N08-N13
      6 1644 1344 (N29W41 L= 54)   376    105        0 (1604/1308-1678/1382) -38.7 W37-W44 N26-N30
      7  729 1233 (N09E21 L=351)  1040    239        1 ( 680/1125- 775/1326)  81.3 E23-E19 N02-N15
      8  939 1208 (N10E08 L=  5)   230     86        0 ( 911/1165- 966/1241) -45.0 E09-E05 N08-N12
      9  416 1122 (S02E41 L=331)   478     85        0 ( 385/1091- 456/1159) -15.9 E44-E37 S03-N00
     10  944  668 (S23E00 L= 12)  2238    647        3 ( 692/ 569-1213/ 768)  32.0 E17-W15 S32-S12
     11  317  612 (S36E53 L=319)  1069    209        0 ( 276/ 527- 380/ 696) -65.6 E57-E48 S41-S30
     12  457  637 (S32E37 L=336)   475    110        0 ( 397/ 624- 512/ 655)   0.0 E42-E31 S32-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