Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1348 1576 (N32W37 L= 30)   515    103        1 (1292/1556-1398/1592) -23.2 W32-W40 N30-N33
      2  703 1458 (N40E13 L=340)   592    123        0 ( 646/1416- 759/1497)  16.7 E17-E07 N37-N41
      3 1066 1252 (N19W07 L=  0)   594    185        3 ( 981/1225-1168/1278) -24.4 W02-W13 N16-N21
      4 1351 1206 (N10W23 L= 16)   564    136        2 (1275/1180-1411/1225) -12.5 W18-W27 N08-N12
      5  562 1008 (N15E29 L=324)   367    146        1 ( 544/ 931- 587/1071)  69.4 E30-E26 N09-N18
      6  544  643 (S07E38 L=315)   244     40        0 ( 530/ 618- 564/ 665)  71.6 E39-E36 S08-S05
      7 1033  584 (S19E09 L=344)   381    145        1 ( 998/ 512-1066/ 640)  45.0 E12-E05 S23-S16
      8  515  540 (S13E44 L=309)   353     67        0 ( 499/ 507- 532/ 578)  78.7 E45-E42 S15-S10
      9  874  316 (S35E29 L=324)   442     58        1 ( 855/ 287- 899/ 340)  45.0 E31-E26 S37-S32
     10 1256  251 (S52E03 L=349)   540    103        0 (1205/ 239-1307/ 265) -18.4 E08-W01 S54-S49
  

    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