Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1243 1670 (N39W14 L=165)   729    204        1 (1155/1639-1352/1693)   0.0 W07-W22 N36-N40
      2  370 1619 (N34E54 L= 96)  2410    344        4 ( 323/1565- 431/1676) -23.6 E63-E47 N29-N39
      3  261 1342 (N15E55 L= 95)  1524    282        7 ( 218/1246- 316/1442) -77.0 E61-E50 N08-N21
      4  741 1265 (N10E18 L=132)   211    116        0 ( 710/1220- 772/1320) -50.2 E20-E15 N07-N13
      5  639 1246 (N09E24 L=126)   210     73        0 ( 594/1235- 677/1256)   0.0 E27-E21 N07-N09
      6  619 1088 (S01E25 L=125)   325     94        2 ( 578/1038- 658/1130) -76.0 E27-E22 S03-N01
      7  336 1031 (S04E46 L=105)   270     41        1 ( 312/1008- 354/1053) -33.7 E47-E44 S05-S02
      8  271 1023 (S05E51 L= 99)   321     49        1 ( 248/1003- 293/1047) -26.6 E53-E49 S05-S03
      9 1598  811 (S14W38 L=188)   785    240        4 (1520/ 700-1663/ 901)  61.7 W34-W42 S21-S08
     10 1398  856 (S13W24 L=174)   437     74        0 (1356/ 827-1444/ 887)  38.7 W21-W26 S14-S10
     11  280  701 (S24E54 L= 96)  1049    248        0 ( 228/ 643- 336/ 745)  24.8 E60-E47 S27-S21
     12  593  570 (S33E30 L=120)   665    188        2 ( 517/ 539- 664/ 612)  13.0 E36-E23 S34-S29
  

    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