Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1315 1681 (N35W39 L=134)  1084    172        1 (1241/1656-1404/1706) -22.6 W32-W45 N31-N37
      2 1572 1602 (N23W54 L=149)   459     78        0 (1535/1575-1611/1623)   0.0 W49-W58 N22-N24
      3 1125 1345 (N19W13 L=107)  1283    409        1 ( 941/1284-1279/1453)   2.4 W01-W25 N16-N21
      4  586 1393 (N34E21 L= 73)  1483    344        2 ( 460/1303- 720/1448) -30.1 E31-E12 N25-N38
      5  173 1361 (N40E71 L= 23)   627     76        0 ( 151/1329- 204/1383)  23.2 E77-E64 N38-N42
      6  156 1269 (N35E69 L= 25)   881     66        3 ( 133/1233- 185/1301) -18.4 E76-E62 N32-N36
      7 1847 1191 (S06W61 L=155)   547    108        2 (1831/1149-1864/1237)  33.7 W58-W64 S08-S03
      8 1747  892 (S21W45 L=140)  1854    370        5 (1677/ 761-1831/1049)  37.4 W38-W55 S28-S13
      9  879  713 (S12E15 L= 79)   514    118        2 ( 862/ 647- 899/ 782)  60.3 E17-E12 S15-S08
     10 1413  686 (S26W17 L=111)   209     40        1 (1395/ 660-1433/ 706) -90.0 W15-W18 S27-S23
     11 1342  525 (S35W09 L=104)   236     52        2 (1316/ 504-1369/ 542) -14.0 W07-W11 S36-S32
     12  349  449 (S18E66 L= 28)   523     83        1 ( 318/ 417- 377/ 491)  80.5 E67-E63 S20-S14
  

    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