Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  787 1862 (N47E49 L=340)  1123    178        2 ( 756/1812- 839/1894)  -5.2 E54-E42 N42-N50
      2  836 1759 (N40E36 L=353)   392    118        0 ( 809/1730- 870/1792) -12.5 E40-E31 N38-N42
      3  550 1482 (N15E41 L=348)  2151    569        9 ( 432/1346- 679/1614) -20.8 E57-E28 N08-N21
      4 1168 1335 (N19W02 L= 31)   560    184        2 (1128/1272-1214/1436) -48.8 E02-W05 N15-N24
      5 1537 1318 (N27W27 L= 56)   644    184        0 (1478/1241-1589/1405) -38.7 W21-W31 N23-N31
      6 1199 1147 (N09W08 L= 36)   319     63        1 (1178/1125-1227/1182) -45.0 W06-W09 N06-N10
      7  544 1155 (S05E32 L=357)   510    142        2 ( 480/1124- 608/1180)  35.5 E35-E28 S07-S01
      8  800  977 (S10E12 L= 17)   566    232        2 ( 690/ 943- 912/1009)  13.0 E18-E05 S11-S07
      9  212  834 (S31E56 L=333)  1363    363        3 ( 148/ 758- 268/ 926)  45.0 E63-E50 S36-S23
     10 1266  550 (S25W27 L= 56)   242     74        0 (1229/ 537-1304/ 563)  11.3 W24-W29 S25-S23
     11  847  527 (S37W00 L= 29)   513    161        2 ( 769/ 500- 927/ 558)  33.7 E05-W05 S39-S32
  

    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