Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1704 1607 (N22W73 L= 18)   904     71        2 (1686/1571-1727/1647)  16.7 W68-W81 N20-N23
      2  793 1421 (N36E06 L=298)   880    245        0 ( 668/1392- 902/1443) -20.6 E15-W01 N31-N38
      3 1133 1378 (N25W15 L=319)   612    236        0 (1027/1311-1211/1409) -51.3 W09-W18 N19-N29
      4 1108 1208 (N16W09 L=313)   449    186        0 (1072/1150-1162/1247)  71.6 W06-W11 N12-N18
      5 1840 1043 (S12W57 L=  1)   558     74        2 (1822/ 999-1858/1081)  53.1 W55-W60 S14-S09
      6 1683 1041 (S08W42 L=346)   284     60        0 (1667/1009-1698/1071)  33.7 W40-W43 S08-S05
      7  880 1041 (N11E08 L=296)   239    116        0 ( 832/1005- 917/1070) -51.3 E10-E06 N08-N13
      8  207  930 (N15E62 L=242)  2141    312        5 ( 110/ 832- 284/1060) -47.0 E78-E53 N08-N23
      9  478 1029 (N17E35 L=269)   270     63        1 ( 450/1007- 509/1050) -26.6 E36-E32 N16-N18
     10  478  890 (N09E37 L=267)   394     82        1 ( 459/ 854- 509/ 929) -78.7 E38-E34 N06-N11
     11  210  702 (N01E68 L=236)   681     88        1 ( 178/ 677- 255/ 729)  -6.3 E75-E63 N00-N02
     12 1627  702 (S27W33 L=337)   265     29        1 (1606/ 679-1646/ 722)  71.6 W31-W34 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