Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1198 1681 (N43E07 L=127)  1103    201        2 (1154/1601-1254/1755) -34.7 E13-W00 N38-N47
      2  744 1712 (N32E40 L= 94)   427     90        1 ( 725/1671- 772/1743) -71.6 E41-E37 N28-N34
      3  794 1494 (N19E26 L=108)   743    225        2 ( 745/1399- 853/1588) -81.9 E28-E23 N12-N26
      4  501 1478 (N11E45 L= 89)   834    201        1 ( 465/1382- 547/1539)  90.0 E48-E43 N05-N15
      5  573 1503 (N14E41 L= 93)   284     63        1 ( 557/1472- 600/1523)  71.6 E42-E39 N12-N15
      6  939 1387 (N16E14 L=120)   299    103        0 ( 898/1347- 973/1430)  74.1 E15-E12 N12-N19
      7  165 1250 (S08E71 L= 63)  1521    271        5 ( 123/1214- 202/1289)  26.6 E78-E65 S10-S05
      8  962  851 (S14W00 L=135)   199     44        1 ( 934/ 829- 983/ 868) -18.4 E01-W02 S15-S13
      9  629  646 (S35E17 L=117)   294     42        1 ( 612/ 620- 655/ 678) -26.6 E18-E14 S36-S33
     10 1250  317 (S40W40 L=174)  1354    288        1 (1144/ 288-1389/ 348)  11.3 W30-W50 S42-S36
  

    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