Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  573 1665 (N36E38 L= 86)  3029    576        4 ( 421/1597- 739/1744) -16.4 E57-E23 N31-N43
      2  536 1418 (N19E34 L= 90)   481    143        0 ( 496/1363- 576/1470) -40.6 E37-E30 N15-N22
      3  472 1404 (N18E38 L= 86)   307     44        1 ( 455/1372- 491/1432)  76.0 E39-E37 N15-N19
      4  396 1306 (N11E43 L= 81)   849    182        3 ( 376/1227- 415/1395) -74.7 E44-E40 N06-N17
      5 1188 1232 (N09W09 L=133)   487    146        0 (1142/1163-1235/1308) -58.0 W06-W12 N04-N13
      6  721 1050 (S04E18 L=106)   253     71        3 ( 703/1019- 742/1086)  76.0 E19-E17 S05-S01
      7 1835  725 (S17W62 L=186)  1363    274        1 (1748/ 630-1903/ 846)  51.3 W55-W68 S23-S08
      8 1757  789 (S14W52 L=176)   984    161        2 (1701/ 755-1816/ 824)  29.1 W47-W57 S16-S11
      9  496  615 (S31E36 L= 88)  1209    380        2 ( 419/ 486- 611/ 732)  41.8 E45-E26 S39-S22
     10 1217  693 (S23W14 L=138)   334     93        1 (1164/ 666-1272/ 724)  29.7 W10-W17 S25-S21
     11  945  567 (S33E04 L=121)   468    109        3 ( 897/ 540-1007/ 591)  26.6 E06-W00 S34-S31
     12  599  325 (S52E38 L= 86)   338     77        1 ( 584/ 300- 620/ 350)  33.7 E40-E34 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