Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  796 1488 (N39E10 L=127)   554    115        1 ( 748/1461- 838/1516)   0.0 E13-E06 N36-N40
      2  245 1264 (N29E60 L= 77)   554     62        3 ( 226/1231- 264/1299) -68.2 E62-E57 N26-N31
      3  937 1211 (N19E03 L=134)   430    131        1 ( 879/1163- 969/1246) -50.2 E06-W00 N15-N21
      4  416 1201 (N24E41 L= 96)   349     57        0 ( 404/1166- 427/1244) -68.2 E41-E39 N21-N26
      5  846 1174 (N18E09 L=128)   292    125        0 ( 829/1113- 871/1231) -90.0 E09-E08 N13-N21
      6 1870  936 (S14W63 L=199)  1575    209        5 (1818/ 838-1909/1021)  34.7 W56-W69 S19-S09
      7  690  943 (N05E22 L=115)   328    133        0 ( 655/ 890- 726/ 998)  45.0 E24-E18 N02-N08
      8  608  521 (S20E36 L=101)   359     59        1 ( 585/ 490- 631/ 544)  21.8 E38-E33 S22-S18
      9  385  470 (S22E60 L= 77)   576    107        2 ( 348/ 431- 415/ 508)  20.6 E69-E53 S25-S19
     10  503  381 (S30E53 L= 84)   975    138        1 ( 446/ 332- 560/ 426)  17.7 E66-E44 S33-S26
     11  924  343 (S38E19 L=117)   478    108        2 ( 864/ 324- 972/ 368)   6.3 E24-E15 S38-S36
  

    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