Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1245 1598 (N31W33 L=165)   242     81        0 (1215/1568-1268/1631)  14.0 W29-W37 N29-N32
      2  957 1533 (N37W12 L=144)   967    280        1 ( 819/1502-1085/1552) -27.9 W02-W19 N31-N40
      3  339 1046 (N24E44 L= 87)   267    104        1 ( 279/1020- 375/1066) -18.4 E50-E41 N21-N25
      4 1569  902 (S16W29 L=160)  1009    355        0 (1446/ 805-1678/1037)  15.9 W18-W39 S18-S12
      5 1201  883 (S07W06 L=138)   928    250        1 (1141/ 766-1253/ 987)  85.9 W04-W07 S14-N00
      6  249  894 (N17E56 L= 76)   546     75        2 ( 215/ 863- 275/ 933) -90.0 E60-E53 N14-N18
      7 1597  806 (S23W29 L=161)   312    118        0 (1564/ 757-1633/ 840)  18.4 W25-W31 S25-S21
      8 1610  716 (S29W28 L=160)   376     80        1 (1582/ 689-1639/ 755) -63.4 W26-W30 S31-S25
      9 1410  582 (S31W11 L=143)   287     44        2 (1384/ 561-1439/ 604) -14.0 W09-W13 S32-S29
     10 1284  473 (S34E00 L=131)  1099    211        2 (1232/ 389-1325/ 556)  57.5 E06-W03 S39-S27
     11  579  453 (S16E45 L= 87)   328     67        0 ( 564/ 417- 593/ 488)  59.0 E46-E42 S18-S13
  

    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