Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  573 1735 (N24E54 L=303)   629    171        1 ( 528/1707- 618/1762)  -4.8 E60-E48 N22-N25
      2  195 1303 (S10E67 L=290)   549     48        1 ( 174/1278- 217/1329)  45.0 E69-E65 S11-S08
      3  291 1203 (S14E52 L=305)   745    143        1 ( 270/1131- 321/1273) -83.7 E54-E50 S18-S09
      4 1169 1196 (N07W03 L=  1)   305    117        0 (1129/1150-1193/1257)  90.0 W02-W04 N03-N11
      5  755 1151 (S06E18 L=339)   256     61        0 ( 724/1122- 791/1174)  53.1 E19-E16 S08-S04
      6 1802 1075 (N21W48 L= 46)  1073    232        1 (1750/ 967-1851/1155) -28.3 W42-W56 N16-N23
      7  275 1055 (S23E50 L=307)   942    225        2 ( 208/1026- 353/1092)   0.0 E57-E42 S24-S21
      8  912  750 (S24W02 L=359)   335     78        0 ( 868/ 730- 953/ 773)  59.0 W00-W03 S26-S21
      9  892  623 (S32W04 L=  2)   680    188        2 ( 807/ 588- 985/ 682)   0.0 E02-W11 S34-S30
     10 1149  517 (S30W24 L= 22)   371     78        1 (1119/ 496-1189/ 553)  -9.5 W21-W27 S31-S28
     11 1702  475 (S14W66 L= 64)   872     98        2 (1666/ 438-1740/ 526)  81.9 W63-W69 S16-S09
     12 1291  314 (S37W45 L= 43)   586    132        1 (1230/ 298-1360/ 331)  45.0 W41-W49 S40-S33
  

    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