Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  513 1556 (N34E43 L= 10)  1306    354        5 ( 429/1456- 578/1622) -34.5 E52-E34 N27-N38
      2 1053 1320 (N21W01 L= 54)  1265    303        1 ( 912/1235-1187/1420) -29.1 E09-W09 N16-N27
      3  210 1351 (N18E67 L=346)  1061    164        1 ( 188/1311- 238/1399) -32.0 E70-E61 N15-N20
      4  707 1295 (N18E22 L= 31)   462    150        2 ( 666/1226- 750/1365) -56.3 E25-E18 N13-N22
      5  567 1250 (N14E31 L= 22)   557    152        1 ( 549/1182- 585/1333) -83.7 E33-E29 N10-N19
      6  337 1184 (N09E49 L=  4)   765    124        2 ( 319/1129- 361/1255) -69.4 E51-E46 N05-N13
      7  386 1062 (N02E43 L= 10)   609    114        2 ( 351/1026- 431/1098) -23.2 E46-E39 N00-N04
      8  282  995 (S03E52 L=  1)   448    170        0 ( 237/ 978- 342/1015)  10.3 E57-E46 S04-S01
      9 1085  265 (S51W11 L= 64)   573    133        0 (1025/ 226-1130/ 296) -26.6 W04-W16 S55-S48
  

    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